Axion as a Cold Dark Matter candidate 
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Here we generally prove that the axion as a coherently oscillating scalar field acts as a cold dark 
matter in nearly all cosmologically relevant scales. The proof is made in the linear perturbation 
order. Compared with our previous proof based on solutions, here we compare the equations in the 
axion with the ones in the cold dark matter, thus expanding the valid range of the proof. Deviation 
from purely pressureless medium appears in very small scale where axion reveals a peculiar equation 
of state. Our analysis is made in the presence of the cosmological constant, and our conclusions are 
valid in the presence of other fluid and fleld components. 

PACS numbers; 14.80.Mz, 95.35.+d, 98.80.Jk, 98.80.-k 



Cold dark matter (CDM), despite its unknown nature, 
became an essential ingredient in the current cosmolog- 
ical studies concerning the large-scale structure forma- 
tion. CDM in a cosmological constant dominated world 
model (often termed the ACDM model) is currently the 
most successful candidate of cosmological models. The 
nature of CDM and the nature of cosmological constant 
(or some other dynamical dark energy), however, still 
remain as fundamental mysteries of present day physi- 
cal cosmology. From early days of dark matter studies 
axion as a coherently oscillating scalar field is widely ac- 
cepted as a candidate for the CDM Confirmation of 
the case using the relativistic linear perturbation anal- 
ysis was made previously ^3, i3i]i for a recent study see 
[jj. In this work we present a more general proof that 
the axion can be regarded as the CDM to the linear or- 
der perturbation. We also derive an effective equation of 
state of the axion which could be important in the solar 
system scale if the system is in the linear regime. We set 
c=l = h. 

We consider the axion as a minimally coupled scalar 
field with V = im^0^. The relevant current scales we 
are concerned correspond to 
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where we set currently H = lOOhkm/ (secMpc). In the 
following analyses we strictly ignore H/m higher order 
terms. But we do not impose any limit on the wavenum- 
ber A:' this is more general than our previous work in Q 
and [J]. We take a spatially flat background with the 
cosmological constant A and the axion; inclusion of A is 
also more general than previous studies. 

We consider the temporal average of the oscillating 
scalar field contributes to the background fluid quanti- 



ties, thus 
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3H<p + m^(f> = 0, 



(3) 



where the angular bracket indicates averaging over time 
scale of order m^^, see Eq. (5) in [3]. Under an ansatz 

(j){t) = {t) sin {mt) + (f)- (t) cos (mt) , (4) 

ignoring H/m higher order terms, Eq. ^ leads to an 
approximate solution 

(t>{t) = a^^/^ [0+0 sin {mt) + 0_o cos (mt)] , (5) 
where (/)+o and 0_o sue constant coefficients. Thus, 
1 



-m^a-^ (Si 



o + 0-o), P = 0, 



(6) 



and the background medium evolves exactly same as a 
pressureless ideal fluid @. 

We consider only the scalar-type perturbations. Our 
conventions are 



ds^ = - (1 + 2a) dt^ - 2al3^adtdx°' 

W [(1 -f 2ip) Sap + 27,„/5] dx^'dx'^, 
1 



(7) 
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= ip + 5p) 5| + -J (^V" ~ -ASf, j a. (8) 

To the linear order, both the vector- type (rotation) and 
the tensor-type (gravitational waves) perturbation equa- 
tions are not directly affected by the presence of the min- 
imally coupled scalar fleld including the axion [bI]. The 
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basic perturbation equations we need are the Raychaud- 
hury equation, the energy-conservation equationj_and the 
momentum conservation equation, respectively 



k + 2Hk + [3H ~ — ]a = A-kG {Sfi + 3Sp) , (9) 



Sfi + 3H {5fi + 6p) = {^i+p)[K- 3Ha - -v ] (10) 

--2^1 -(11) 
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The temporal average of oscillating scalar field con- 
tributes to perturbed fluid quantities as 0] 



= {(pSc/) — (pa + m' 
Sp — {(l)S(j) — (jy^a — m' 



-[ii + p)v 



cr = 0. 



(12) 



For the scalar field we have the equation of motion [6] 

S(j) + 3HS(}) + -^5(j) + Va,a,5(I) 

= <P{K + a) + (24> + 3H<j>^a. (13) 

The above set of equations is spatially gauge-invariant. 
But we have not taken the temporal gauge (hypersurface) 
condition which will be used as an advantage in handling 
the mathematical analyses. 
For S(j) we take an ansatz 

6(j){k,t) ^ 5(f)+{k,t) sin {mt) + 5 (j)^{k,t) cos (mt). (14) 

Using Eqs. ^ and HI]), Eq. gives to the leading 
order in H/m 
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TO (0+0(50- - 0-Q(50+) . (15) 



We have not taken the temporal gauge condition yet. 

The comoving gauge takes v = as the temporal gauge 
(hypersurface) condition; in the presence of additional 
fluid, like dust, radiation, neutrinos, etc., our gauge cor- 
responds to the axion-comoving gauge. Equation (lll|) 
gives 



Sp 



(16) 



Equation (fTUl) . and Eqs. © and (|lip . respectively, give 
S = K, (17) 



k + 2Hk = AttGhS 



k"^ Sp 



(18) 



where S = Sfi/ ^. Combining these equations we have 



S + 2HS - AnGfiS - 



fc^ Sp 



0. 



(19) 



This equation is well-known in the Newtonian context, 
see Eq. (15.9.23) in [71], and also valid in Einstein grav- 
ity context in the limit of negligible pressure [§] . In the 
case of CDM we effectively have Sp = 0. However, the 
equation of state relating Sp with Sfi in the case of axion 
is not determined at this point. The equation of state of 
axion will follow from the perturbed equation of motion. 
Under the comoving gauge, Eq. gives 



(20) 



Relation between Sp and Sfi can be determined through 
a using Eq. (fT3|) . We solve Eq. (fT3|) strictly to leading 
order in H/m as the solution for the background is valid 
to such an order. Using Eqs. (fTH]) and (fT7| . Eq. (fTS]) gives 
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Equation gives 
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0+0 V 4 rn^a'^ 
Therefore, we have an equation of state 
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The effective sound speed becomes 




-1/2 



(21) 



(22) 



(23) 



(24) 



which shows an interesting scale dependence; for 
k/(rna) <C 1 the time dependence Cs oc 1/a is the same 
as in the ordinary matter dominated medium, see later. 
Equation (fT9|) becomes 



S + 2HS 



^AttGh 



1 k^ 
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4 TO^a^ 1 
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(5 = 0. (25) 



We note that we only have assumed (H/m)^ ^ 1, but 
have not assumed any condition on [k/aHY'. Thus our 
equations are valid in all scales. Even in the small scales 
where the pressure gradient term has a role, in gen- 
eral we have k^/{m'^a'^H'^) ^ k"^ /{rn^a?), thus ignoring 
k"^ / {m?a?) term we have 



5 + 2H5 - { AttGh 



1 k^ 



4 m?a^ 



0. 



(26) 



This form was presented in Eq. (20) of [i^ based on the 
zero-shear gauge (often termed the conformal Newtonian 
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gauge or the longitudinal gauge). Although Eq. ([26]) . 
more precisely Eq. (PS)) . in our comoving gauge is valid 
in all scales, in the zero-shear gauge it is valid only for 
k/{aH) 3> 1; the general density perturbation equation 
of pressureless medium in the zero-shear gauge is quite 
complicated, see Eq. (27) in Q. The competition be- 
tween gravity and pressure gradient terms in Eq. (I26p 
gives the Jeans scale; assuming a flat axion dominated 
model we have currently 



Xj = — = 2n (l67rG/im2) 



1/4 



= 5.4 X lO^^cmh 



-1/2 ^ 



-1/2 



(27) 



which is quite small corresponding to the solar-system 
size for m ~ 10~^eV, see also On scales larger than 
Aj, thus effectively in all cosmologically relevant scales, 
the axion fluid can be regarded as the pressureless ideal 
fluid. Therefore, the axion is justified as a CDM candi- 
date. 

Equation ([26|) can be analytically solved for A = 0. 
For the background, we have a oc t'^^^, H = 2/(3t), and 
II = l/{6nGt'^). Equation p6)) has an exact solution 

l2 /^2/3\2 

5±^t-'/'j^,M3At-'/% A^ — ^—j . (28) 

This can be compared with an exact solution in mat- 
ter dominated era with an equation of state p (x ^Ti 
0. For 7 = 5/3 we have the solution in Eq. (gS]) 
with A = Cskf^/^ /a, see Eq. (15.9.4f) in [7|]; notice that 
Cs oc f/a even in the matter dominated medium. Thus, 
again we can identify = fc/(2ma) as an effective sound 
speed of the axion fluid. 

For fc^/ [raaH] -C 1 and A = we have a perturbative 
solution 



5{k,t) ^ c+{k)t^'^ 



1 



3k^ (el^ 
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(29) 



This coincides with the solution in Eq. (25) of [1]. Our 
previous proof of the axion as a CDM candidate pre- 
sented in [31 was made based on this and other solutions. 



For a clear presentation, in this work we have consid- 
ered a flat background composed of a single axion compo- 
nent but including the cosmological constant. Our anal- 
ysis can be easily extended to non-flat case as well as in 
the presence of additional fluids and fields. We can show 
that axion behaves as a CDM even in such more realistic 
situations. For example, in the presence of other com- 
ponents the equation of state in Eq. is valid for the 
axion component. 

Equation (j27p shows an axion Jeans scale where linear 
perturbation of axion fluid becomes stabilized (oscillates) 
under that scale. The axion fluid shows peculiar equation 
of state and effective sound speed presented in Eqs. P5|) 
and . Although our neighborhood in the solar system 
is in a significantly nonlinear stage it is curious to see 
the possible observational signature of axion fluid based 
on its contribution to effective pressure. This is left for 
future investigation. 

Recently we have shown that even to the second-order 
perturbation, the density and velocity perturbation equa- 
tions of the zero-pressure medium in Einstein's gravity 
exactly coincide with the ones in Newton's gravity: we 
call this a relativistic/Newtonian correspondence p^ . 
The relativistic/Newtonian correspondence to the lin- 
ear order can be found in Eq. (|19p with vanishing pres- 
sure. In the present work we have shown that the axion 
properly treated in relativistic perturbation theory be- 
haves as a pressureless fluid in cosmologically relevant 
scales, thus justified as a CDM candidate: we may call 
it the axion/CDM correspondence to the linear order. 
Whether such a correspondence continues even in the 
case of second-order perturbation is an interesting open 
issue at the moment. We will address this important 
issue in a future occasion. 
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